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We present a technique to calculate the transport properties through one-dimensional models 
of molecular wires. The calculations include inelastic electron scattering due to electron-lattice 
interaction. The coupling between the electron and the lattice is crucial to determine the transport 
properties in one-dimensional systems subject to Peierls transition since it drives the transition 
itself. The electron-phonon coupling is treated as a quantum coherent process, in the sense that no 
random dephasing due to electron-phonon interactions is introduced in the scattering wave functions. 
We show that charge carrier injection, even in the tunneling regime, induces lattice distortions 
localized around the tunneling electron. The transport in the molecular wire is due to polaron-like 
propagation. We show typical examples of the lattice distortions induced by charge injection into 
the wire. In the tunneling regime, the electron transmission is strongly enhanced in comparison 
with the case of elastic scattering through the undistorted molecular wire. We also show that 
although lattice fluctuations modify the electron transmission through the wire, the modifications 
are qualitatively different from those obtained by the quantum electron-phonon inelastic scattering 
technique. Our results should hold in principle for other one-dimensional atomic-scale wires subject 
to Peierls transitions. 
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I. INTRODUCTION 

Developments in nanofabrication, including both 'bottom-up' approaches and 'top-down' methods are focussing 
renewed attention on the properties of molecular-scale entities and their potential as electronic device componentscl. 
One of the most basic theoretical questions that can be asked in this context is: what determines the conductance of 
a molecule if it is used as a current-carrying element bridging two reservoirs of differing electron chemical potential? 
This question now has an immediate relevance for experiments in which-such conductances are measured, using either 
(i) scanning probe tips for individual malecules adsorbed on surfacescTEl, for molecular wires adsorbed at^step edgesQ, 
embedded in self-assembled monolayerg|3ij,pOi'i-(ii) (macroscopic) electrodes obtained by nanolithographyO Ej or from 
a mechanically controllable break junction£J"E3.p-| 

Since the seminal work of Aviram and Ratnerll3 concerning the electron transfer rate between acceptor and donor 
groups linked by a conjugated molecular bridge, numerous theoretical studies on electron transfer and transport 
through molecular systems have been performed. In the following, we briefly review some contributions on the 
electron transport through a single organic molecule (or a few molecules) whose ends are connected to electron 
reservoics—Calculations of the electronic tmiismission through such systems have been done for purely one-dimensional 
modelsll3t3 and two-dimensional modelaH3c3. More realistic descriptions of the electrode/molecule system have also 
been developed. Combining elastic electron scattering theories with three-dimensional tisht-iiiiLding-like Hamiltonians, 
models,have been developed for molecular wires connected to two semi-infinite surfacesEjoOLJ or to two semi-infinite 
'rods'EJ, or to cluster-like leads where imaginary parts pam^pjtroduced in the Hamiltonian to take into account the 
fact that electrons can leak into the metallic reservoirsEEloEj. Within a framework equivalent tOj-the latter model, 
calculations have been extended to the Hartree-Fock level for a molecule attached to gold clustersEfl. More recently, 
density functional theory has been applied to a molecular wire (described by atomic pseudopotentials) connected to 
two jellium surfacesL3. 

These theoretical studies have clarified the importance of three major points crucial for the transport properties 
of the molecular wire connected to the reservoirs. First, they have shown the importance of the electronic and 
chemical interaction between the ends of the molecular wire and the reservoirs. The larger the Hamiltonian matrix 
elements between the delocalized electronic states of the electron reservoirs and those molecular electronic states 
that extend along the wire, the better the conductance properties will be; furthermore, these matrix elements should 
be large compared witii-|the characteristic intra- molecular Coulomb interaction between electrons in order to avoid 
the Coulomb blockadeE£l. Second, the theories show that in the limit of small applied voltage and away from the 
Coulomb blockade regime, the transport is dominated by charge carrier tunneling inside the HOMO-LUMO gap of 
the molecule. This gap is another crucial parameter for control of the conductance of the wire. The smaller the gap is. 
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the larger the tunneling transmission will be. More generally, the gap of a molecule depends on the chemical nature 
and atomic structure of the system. This gap can also be modified by the electron-electron interactions or by a change 
of the structure of the molecule due to (z) external forces, (ii) (thermal) lattice fluctuations or {Hi) electron-lattice 
interaction. Third, the calculations highlight the importance of the position of the molecular electronic levels with 
respect to the Fermi levels of the reservoirs in the presence of an applied voltage^nand the related issue of wJaere 
the pptcntial drops occur inside the junction. This has been done both empiricallyE3 and by using approximatecj or 
exactE3 self-consistent schemes. 

Although the above theoretical work has shed light on several important physical processes for the transport in 
molecular wires and has matched to a certain level of accuracy with some experimental measurements, all the models 
previously cited are based on elastic electron scattering through the rigid lattice of the wire. However, in such highly 
confined electron systems, the coupling between electron and other excitations (phonons for instance) is strongly 
enhanced because of the size of the system and its quasi one-dimensionality. This makes the rigid lattice approximation 
questionable^particularly so since a one-dimensional metal is generically unstable to a Peierls transition at low 
temperaturet3. In an infinite system, such a transition typically produces a semiconductor in which the states near 
the band extrema are very strongly coupled to distortions of the system; in a conjugated organic molecule, the 
corresponding phenomenon is a strong coupling of the 7r-electrons occupying the HOMO and LUMO states to the 
bond-alternation pattern. This coupling means that the low-lying states of a charged molecule (via which any net 
transport of charge through the molecule must proceed) involve-apJntimate coupling of electronic and lattice degrees 
of freedom, to produce excitations such as polarons or solitonsESLJ. These coupled excitations can be thought of as 
conspiring to lower the energy gap locally around a charge carrier when it is introduced into the system. Such uolaxcpjc 
and solitonic phenomena have been studied in bulk or thin film samples of conducting polymers for decadesLJtj CJ. 

The importance of this electron-lattice coupling means that the conventional manner of introducing lattice vibrations 
within a Landauer-type approach to conductance, as an extrja-bijoadening of the electronic levels (extra imaginary 
part in the corresponding Hamiltonian, see for example RefsHj'L3), is not sufficient to describe the coherent lattice 
distortion due to charge injection. To our knowledge, the explicit natuce of the distortion accompanying charge 
injection has only so far been partially addressed in two simplified limitsEl. In the first case, a molecular wire was 
treated as a rigid lattice in which a static soliton-like defect is prcscntUJ. Although it possesses a mid-gap electronic 
state, this model does not permit the study of the dynamics of formation and transport of charge-induced lattice 
distortions. In the second case, tJie atoms of a conjugated molecule were assumed to respond classically to the 
injection of a electron wave-packetCj, via forces calculated from expectation values of the electron wave-packet and 
other electronic states. In this model, the lattice is able to respond to the injected charge, but not in the physically 
correct manner: within a wave-packet approach to tunneling, only a small part of the electronic charge enters the 
tunnel barrier. Therefore the lattice responds with probability unity to a small fraction of the charge of the injected 
particle, rather than responding with a small probability to the total charge of the injected particle. 

The only way to overcome these limitations is to perform transport calculations in which the full dynamical cor- 
relation between charge carriers and quantum phonons is retained. We report the results of such calculations in this 
paper. In order to focus on this particular mechanism for charge transport, we use a jSipiple tight-binding model 
of a conducting polymer (the Su-Schrieffer-Heeger (SSH) model for trans-polyacetyleneES) that does not explicitly 
include any electron-electron interactions. However, our calculations cover a range of transport regimes that includes 
tunneling transport (virtual electrons) and resonant transport (where there is sufficient energy to inject real electrons 
into the system). In contrast to the more usual 'phase-breaking' approaches to the electron-phonon interaction in 
transport problems, we explicitly retain the phase coherence between elastic and inelastic processes. In this paper we 
concentrate on systems where the boundary conditions on the molecular wires force them to be semiconducting in the 
zero-bias limit. Within the SSH model, this correponds to molecular chains containing an even number of monomers. 
Calculations on odd-length chains, which incorporate mobile solitonic defects with associated mid-gap states that 
make an additional contribution to the transport, will be reported separately. 

The paper is organized as follows. In Section ||, we present the multichannel scattering technique used to calculate 
the transport properties of the molecular wires. This section also includes a detailed analysis of different approxima- 
tions used for modelling the molecular wires (involving harmonic phonons) and for reducing the computational cost of 
the calculations (involving reduction of the parameter space). The results obtained, in the limit of low temperatures, 
for the response of the molecular wires to charge carrier injection are given in Section HI. We show how the lattice is 
distorted by the injection of a tunneling electron and how the coherent coupling between the tunneling electron and 
the quantum phonons affects the transmission properties through the wires. We also compare our transport results 
with the effect of straightforward static fluctuations in the harmonic lattice, excluding the dynamical correlation of 
the electrons and the phonons. Finally we summarize the most important results and propose further developments 
of the present study (Section |^). Additionally, in an appendix we recall briefly the methods used to get the ground 
state and the harmonic phonon modes from the original SSH model. We also derive in the appendix the quantuHi 
electron-phonon Hamiltonian used for the molecular wires. A brief account of part of this work has already appearedcZl. 
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II. MODEL 



We are interested in modelling the coherent electron (or hole) transport through a finite size system (the molecular 
wire) connected to two leads which inject or collect the charge carriers. Within the wire, the charge carriers interact 
with the atomic motion which originally drives the Peierls transition in the molecule. 

The interaction between the ends of the molecular wire and the leads is supposed to be strong enough to permit a 
good overlap between the electronic states of the wire and the surface electronic wave functions of the leads. In most 
of the practical applications, molecular wires end in "active" chemical groups, like Ihiol (S-H) for example, which are 
known to react easily in presence of a gold surface to form chemical Au-S bondsc3. We therefore assume that the 
electron transfer rate at the molecule/lead interface is such that we can consider the electron (hole) transport as being 
a coherent process throughout the nanojunction, rather than a sequential incoherent two-step process. 

In the coherent transport regime, a stationary state wave function scattering technique can be used to calculate 
the electron transfer through the molecular wire. Since we assume that the basis sets used to describe the leads and 
the molecular wire form a complete set, there are basically two ways to solve tlie scattering problem for a single 
incident charge carrier. The technique is reminiscent of the Lowdin transformationca. If one projects out the basis set 
associated to the molecular wire, the problem is reformulated as a "single impurity" wit h . pii -site energies and coupling 
matrix elements to the leads depending on the injection energy of the charge carrieill3H. If one choose to project 
out the basis set associated to the leads, one can effectively remove the leads from the problem. This technique is 
identical to the embedding technique where a (finite gize)_effective Hamiltonian describing the region of interest is 
obtained by introducing complex embedding potentialsLElo. The embedding potentials characterise the matching of 
electronic spectrum of the wire to the continuum of states of the semi-infinite leads. These potentials also depend on 
the charge injection energy. 

As we wish to obtain the response of the molecular wire to charge injection as well as the transport properties 
through the junction, we choose the embedding approach to solve the electron (hole) transport in the system. For 
this, we use a technique whiejMpermits us to map the many-body electron/phonon problem onto a single-particle 
problem with many channelsEj'E^. 

In the remain of this section, we present the basis of the many-channel scattering technique and discuss the 
different approximations introduced to reduce the computational cost of the calculations. Details of the construction 
of thC|4iuaiitum Hamiltonian for the molecular wire from the model originally proposed by Su, Schrieffer and Heeger 
(SSII)l3'L^ are given in Appendix 

Starting from this model, we have derived a quantum electron-phonon Hamiltonian: 

^^en elfin + ^ TT-L^q a\aq + ^ lqnm{a\ + Og) cJ^C^ , (1) 
n q g,n,m 

where c\ creates an electron in the nth one-electron state of the molecular wire with energy e„ and aj creates an 
excitation in the gth eigenmode of vibration (phonon) of the molecule with energy hujq. The Hamiltonian Eq. (|l|) 
goes beyond the Holstein and Frohlich model for the electron-phonon (e-ph) interaction, in the sense that the electron 
couples to different non-local eigenmodes of vibration, each mode having a different frequency. The electron-phonon 
coupling is linear in the phonon field displacement and involves electronic transitions via a general form for the e-ph 
coupling matrix elements "fqnm- 

The electronic eigenstates and eigenvalues are determined self-consistently with the atomic configuration for the 
ground state of the neutral dimerized molecular chain taken to be the reference system (App. From the atomic and 
electronic structures of the reference system, the phonon modes and frequencies are calculated within the harmonic 
approximation (App. The e-ph matrix elements are derived from the SSH model by expanding the atomic 

displacements induced by adding a charge onto the phonon modes of the neutral molecule (App. H). We have 



checked the accurary and the validity of the harmonic approximation for the phonons (see below section [IB) 



A. Multichannel scattering technique 

We are interested in solving the problem of electron transport through a 'nanojunction' within a two-terminal 
device. We are mostly interested in the coherent regime for the electron transport, that is the regime where no 
random phase breaking is arbitrarily introduced between different electron scattering states. Furthermore, we wish to 
use a formalism that can treat (i) different transport regimes (pure tunneling, resonant tunneling, eventually ballistic 
transport) on an equal footing or in a transparent way, and (ii) the coupling of an electron with other degrees of 
freedom within the 'nanojunction'. In principle, to study the eleotsopic transport in such open systems, one would 
have to deal with the non-equilibrium Green's functions formalisn£3~E3. 
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In this paper, we use a model which maps a many-body ocniJem (to be accurate a one electron/many bosons 
problem) onto a single-particle problem with many channelSo. In a such model, one deals directly with the 
multichannel scattering states, although the formalism can be reformulated in terms of Green's functions. This 
multichannel scattering technique has already been used to study electron transport through onc-dimcnsi(j )D|£il . p pdcls 
of (i) double-barrier resonant tunneling junctions with electron coupled to localised single phonon modcESEJH, (ii) 
Holstein phonon model injpwjsence of an electric fieldEj, (iii) mesoscopic structuresEj and Aharonov-Bohm rings with 
on site phonon couplingL^E^, (iv) tunneling barriers with the electron coupled to surface plasmon modesEj. More 
recently such a techruque has also been used to study inelastic electron tunneling through small molecules in a STM 
tunneling barrierE3E3. 

We start with the following heterojunction: a molecular wire containing Na atomic sites, described by the Hamil- 
tonian in Eq. (|l|) is connected to ideal one-dimensional right (R) and left (L) metallic leads, whose Hamiltonians 
are 



+00 



Hr= J2 d\di + (3^{d\di.i + d\_^di) (2) 



and 



via the coupling matrices. 



and 



l=N^ + l 





ffL = V CL d\di + /3L(dJrfi-i + d\_^di) , (3) 



Tk = WR(rf]v +iCw„ + c]v rfjv„+i) (4) 



Tl = wl(4ci+cI'^o) • (5) 

The operators dj (di) create (annihilate) an electron on site / inside the leads, with on-site energy el.r and nearest 
neighbours hopping integrals /3l,r, the operators c| (q) are the corresponding operators for an electron on site i within 
the molecule, and ur and Vh are the hopping matrix elements between the ends of the molecule and the right and left 
leads respectively. Within the molecular wire, the transformation from site representation to eigenstate representation 
is easily performed knowing that cjj = '^f=i Z'^ii) cj, where Zn{i) are the components of the nth electronic eigenstate 
of the wire (see Appendix H). 

The procedure for mapping the problem onto a single-particle system with many channels is performed by writing 
the total scattering wave function |^'(£^)), for the total energy E of the electron-phonon system as 

l*(£^))-EE"Un,}(^)IU^J), (6) 

where the basis set used to expand the scattering waves is defined (in the case of electron transport) as |^, = 
4na^4=|0) (for I < I < Na) and = dlllJ^j=-\0) (for other values of 0, |0) being the vacuum state 



and {rig} being the phonon occupations. The vacuum state |0) is taken to be the neutral ground state of the system, 
with a definite number of electrons in each of the left lead, right lead and molecule. The electronic states we consider 
therefore involve adding a single electron to this neutral state; the added electron may be anywhere in the system 
(in the left lead, the right lead, or the molecule). For hole transport we use an identical basis, except that electron 
creation operators are replaced by annihilation operators. By writing the wave function coefficients as in Eq.(|^), no 
explicit separation between the electronic and phonon degrees of freedom has been assumed. 

As far as the electron (hole) propagation is concerned, each different channel is associated with a different set 
of phonon occupation numbers {ug}. The total wave function \'ii{E)) is the eigenstate of the total Hamiltonian 
H = Hi^ + Ti, + H„ + Tu + Hu, H\'^{E)) = E\'ii{E)), with the full scattering boundary conditions applied {i.e. an 
incident electron or hole with the molecule in a given vibrational state). 

In absence of dissipation, the total energy E of the system is conserved during the scattering process, i.e 

E = Ein+'^Ug hujg = Eont + E ™9 ' 
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where Ei^ is the energy of the incoming electron and {ug} is the initial set of phonon occupation numbers. Eout is 
the energy of the outgoing reflected or transmitted electron with the corresponding set of phonon occupancies {rUq}. 
For inelastic scattering processes {rig} ^ {mq}; for elastic scattering the phonon distribution is conserved. 

The wave function coefficients a; take an asympotic form inside the leads. The form corresponds to propagating 
Bloch waves inside the left and right leads whose amplitudes are the reflection ''{m,} transmission ^jm,} coefficients 
of the electron in the different channels. For an injected electron from the left lead (for example), we have 



a;,{m,} = e'''t"'!}' (5{„j^{„j + r{,„ j e '''^"'!}' (8) 



inside the left lead (/ < —1) and 



ifc^ I 

ai,{m,} = i{m,} e (9) 



inside the right lead (/ > Na + 2). The k^^y are the dimensionless wave- vectors of the Bloch waves in the different 
channels. For a given total energy E, the wave-vectors depend on the phonon occupation numbers. Figure |^ shows a 
simplified sketch of the multichannel technique for different phonon excitations inside the molecular wire. 

The solutions of the Schrorindger equation {l\H\'^{E)) — E{l\'i'{E)) inside the leads give the dispersion relations 
for the electron wave-vectors inside the different channels, Ein = el + 2/3LCOsfc/ i for the incoming wave and 

Eout = El + 2/3l cos k^rng} outgoing reflected wave and to Eout — cr + 2/3r, cos k^^^^ for the transmitted wave 

to the right lead. 

In this paper, we report calculations valid in the limit of low temperatures assuming that fceT <C hujg. Therefore 
we take for the initial set of phonon occupation numbers, the set corresponding to all phonon modes in the ground 
state {riq} — {0}. The dispersion relations become simply E = el^r -I- 2/?l cos fej^^^ for the elastic channels and 

E = cl.r + 2/3l,r cos fc^^ J -I- '^^mqhujq for the inelastic channels. Note that in this case, the total energy E also 
represents the electron injection energy. 

Solving {l\H\^{E)) — E{l\'^{E)) at the interfaces of the heterojunction is the next step to perform in order to 
resolve the unknown wave function coefficients inside the molecular wire. For I — and I — Na + 1, the relations 
between the reflection and transmission coefficients with the wave function coefficients at the ends of the wire are 
obtained: 

= ^ "Wo,{m,} e f'"9> , (10) 

and 

r{mg} = -^{0},{rn,} + ^ > (H) 

for the elastic {0} and inelastic {rriq} ^ {0} channels. 

Finally, solving the Schrodinger equa-tion on sites I G [l,iVa] permits one to effectively remove the leads by intro- 
ducing complex embedding potentialsO. Then the solution of the full scattering problem is obtained by solving the 
following complex linear system 

{E-H^~ Y}^{E) - Y}^{E)\ \a{E)) = \s{E)) , (12) 

where is the molecular wire Hamiltonian Eq. (|l]), the components of |a) are the scattering wave function coefficients 
inside the molecular wire expressed in the molecule eigenstate representation, i.e. a„_{„^} = Zn{i) ai^|„^}, and 
E'"'^ are the embedding potentials due to the left and right lead respectively. The embedding potentials are diagonal 
matrices in the |n, {riq}) basis set with components 

Km^^) = Zn{l) gl^^yiE) «L ^„(1) , (13) 

and 

E^,{„^j(i?) - Z^{Na) vn 5f„,}(£^) ^^R Zn{Na) , (14) 

L R, 

where g^^ ^ are the surface Green's function of the isolated left and right leads for the different channels given by 

= eM^k\f^yiE))/pL..R ■ (15) 
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Finally \s{E)) represents the source term (i.e. the injected electron or hole at energy E) with components given by 



,{n,}iE) = S{o},{n,} (-2i WLsinfc[oj) Z„(l) 



(16) 



In the present case, the boundary condition chosen for the source term corresponds to the injection of the charge, 
into the elastic channel {0}, from the left lead towards the right lead. 

The solution of Eq. can be obtained by several means using algorithms for sparse matrices. In the present 

work, we explicitly separate the real and imaginary parts of vectors and matrices as follows 



(17) 



where 5fte and 5m denote the real and imaginary parts of the different quantitites and E = E'^ + S^ . The matrix 
in the left-hand-side of Eq. ( p^ ) is therefore real and symmetric. We then use stancktrd conjugate gradient (CG) 
technique to solve the linear system A\x) = \b) where A is a real and symmetric matrixEZI. An efficient algorithm has 
been devised for computing the products Hvj\xi) generated during the iterative CG steps. It is based on an optimal 
adressing of the vector components and uses the selection rules for the 7g„m matrix elements (see Appendix A3). 
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B. Isolated molecular wire : harmonic phonons and reduced parameter space 

We have shown that the solution of the scattering problem is obtained by solving Eq. ( [l2| ) for the value of the 
wave functions inside the molecular wire. Assuming a truncated phonon space up to a finite number of excitations 



n 



occ 



, the size of the basis set is given by A'size = Ne x (n™^'^ -I- l)"?*" with N,, (A^ph) being the number of electronic 
states (phonon modes). Even for relatively short wires, the size Ngizc of the basis set quickly becomes too large for 
tractable numerical calculations and/or reasonable computing times. For example for Na = 20 atomic sites (A^ph = 19 
acoustic and optic phonon modes) with only n^^^ — 2, we obtain Agizc ^ 10^. In the following we show how to reduce 
the basis set size and prove the validity of the approximations introduced. The corresponding Hilbert space can be 
reduced by (j) considering only valence (conduction) band electronic states (this will correspond to hole (electron) 
transport respectively), by {ii) considering a limited but sufficient set of phonon modes and finally by (iii) using only 
a few excitations in each phonon modes (truncated harmonic oscillator approximation). 

In order to determine which phonon modes are mainly contributing to the charge-induced deformation of the chain, 
we calculate the ground state atomic configurations for a neutral chain and for a chain charged with one additional 
electron uf. Note that because of the charge-conjugation symmetry (which holds exactly for the SSH model) adding 
an extra electron or removing an electron {i.e. adding a hole) produces exactly the same lattice distortion. Such 
a lattice distortion is known as an electron (or hole) polaron (c./. Fig. |^). The lattice distortion is then projected 
onto the harmpaic eigenmodes of vibration Vq of the neutral chain and the corresponding Huang- Rhys factors Sq are 
determined byO 

Sq = ^ ' , (18) 

where A, = J2i ^gl*) The Huang-Rhys factors give the averaged number of quantum phonons that would be 

needed to achieve the corresponding elastic energy of the lattice distortion. We have checked that the largest values 
for tbcjse factors are obtained for the lowest energy (longest wave-length) optical phonon modes of the molecular 
wiresEj. In particular, the optical phonon mode having the lowest energy also has the most important contribution 
{i.e. Sq > 0.8 as can be seen for different chain lengths in Table |). 

From these results, we can already reduce the number of phonon modes necessary by considering only the longest 
wave-length optical modes to describe the lattice deformation induced by adding a charge into the chain. As typical 
set of these modes is shown on Figure |[ 

Furthermore, we can also ch eck the validity of the harmonic approximation used to determine the phonon modes of 



the wire. Note that in section |A 2\ the elastic energy is expanded up to second order for small displacements around 
the equilibrium atomic positions in order to obtain the dynamical matrix from which the eigenmodes of vibration are 
determined. The deformation of the lattice due to charge addition is not purely harmonic because of the distance 
dependence in the electron hopping matrix elements, cf. Eq. ( |Al[ ). However there are good reasons to believe that the 
potential surface can be fairly well described by the harmonic approximation to the neutral chain. To confirm this, 
we compare the deformation energy AE with the harmonic distortion energy i?harm — ^q ^^I^q'^q — ^q Sq. 
The deformation energy is obtained from AE = Eo[u^] — Eo[u^] where i?o[Ui] is the selfconsistcnt total energy of 
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the SSH Hamiltonian for a neutral chain with the corresponding atomic positions £^o[i''i] is the non-selfconsistent 
total energy of the neutral chain where the atomic positions u1 are taken to be those of the charged chain. In -Eoi^ti] 
the effects of the imposed lattice distortions on the electronic spectrum are taken into account. Typical values of A_E 
and -Eharm are given in Table ^ for different wire lengths. For short chains, the values of AiJ and -Eharm are almost 
indentical (less than k, 5% difference). For longer chains, the difference between IS.E and -Eharm increases but never 
exceed « 15%. We also show below that the harmonic expansion of the elastic energy is sufficient to describe the 
formation of a static polaron by adding a permanent extra charge in the molecular chains. 

Now we turn on the reduction of the Hilbert space in relation to the electron states. We want to check the validity of 
using only half the electronic spectum on the values of the^ atomic displacements induced by adding a charge (electron 
or hole) into the molecular chain. As mentioned in Ref.cZl, we consider the action of the model Hamiltonian Eq. (|l|) 
on the {Na ± l)-electron Hilbert space obtained by adding a charge into the neutral chain. It is thought that it is 
sufficient to consider only the {Na ± l)-electron because for molecular wires strongly (electronically) coupled to the 
leads, the mean time between charge passages is ~ 10~^ s (for a corresponding current of 1 pA), a value orders of 
magnitude bigger than a typical residence time (« 10^^'^ s). Then we project out the addition of a charge into the 
electronic states by working with the electronic eigenstates representation, i.e. the sums over the eigenstates will 
only include the occupied valence band states when adding a hole (and later for hole transport) and only the empty 
conduction band states when adding an electron (for electron transport). In the following, we show that the static 
lattice distortions due to adding a charge are well reproduced by considering only one half of the electronic spectrum. 

We start from the ground state of the reference system (neutral chain with atomic positions u^). The lattice 
distortions due to charging are expanded onto the harmonic phonon modes Vq(i) as Ui = + Vq{i)l!^q. Then the 
influence of the lattice distortions on the electronic Hamiltonian = e„(5„m are taken into account by introducing the 
corresponding e-ph coupling off-diagonal elements 7g„,„ in H'^ . The electronic Hamiltonian of the isolated molecule is 

^nm ~ ^nm ~^ ^nm' ~ ^n^nm + ^ ^ ^q^qnm j (19) 

q 

where the matrix elements "fqnm are given by Eq. ( |A9| ) in Appendix and Aq is the dimensionless displacement 
Aq = Aq^Jfij {2MLUq). The total energy of the distorted lattice is 

i^haif ({A J) = J2 l^H^l + Tr[p^'if<='({A J)] , (20) 

where p'^^ is the electronic density operator. To find the corresponding ground state, i?haif ({A^}) has to be minimised 
versus the classical lattice phonon displacements {Aq}. 

Within the half spectrum approximation, the trace in Eq. ( pO[ ) runs only over the conduction (or valence) band 
eigenstates when we consider the wire charged by an extra electron (or hole). Adding a charge to the system 
involves only the LUMO or HOMO electronic state for an electron or hole respectively, therefore the functional 
to be minimized is actually ii^haif ({A^}) = iMw^A^ + Ae({Aq}) where Ae = {ipe\H'^^{{Aq})\ipf.) is the lowest 
(or highest) energy eigenvalue of the electronic half-space Hamiltonian i/°'({Aq}) and corresponds to the LUMO 
(or HOMO) as modified by the atomic distortion. The minimisation of i?haif({Aq}) is obtained when the forces 
Fq — MujqAq + (fel^qlfe) \/ {2Mujq)/h arc zcro for all the phonon modes q, 7^ being the e-ph coupling matrix with 
components ^qnm- With this procedure, we can also study the contribution of the different phonon modes q and 
check the validity of using only a limited number of optical phonon modes to create the distortions as we have already 
suggested from the values of the Huang-Rhys factors. 

Figure ^ shows the lattice distortions calculated within different approximation from the ground state of a charged 
chain. The lattice distortions are best represented by the staggered difference di between adjacent bond lengths. The 
quantity di = (— l)*(ui+i — 2ui + Ui-i) is known as the dimerization. A constant dimerization pattern indicates a 
perfect bond length alternation in the chain, while a decrease of the dimerization indicates a deformation of the bond 
lengths (i.e. an increase of the short bonds and a decrease of the long bonds). These lattice distortions are localized 
around the charge added and are characteristic of the polaron defect in the molecular chain. The general shape of 
the dimerization pattern in Fig. ^ indicates that a static polaron has formed in the chain and it extents over several 
atomic sites. 

The dimerization patterns obtained from the original SSH model are shown on Figure ^ with the dimerizations 
patterns calculated by considering half the electronic spectrum and classical phonon modes (all the modes and the 
reduced set of long wave-length optical modes). We can see that the lattice distortions representing a static polaron 
defect in the chain are very well reproduced by considering only the long wave-length optical phonons (for example 
those shown on Fig. ^ for the Na = 100). The discrepancies between the dimerization amplitudes obtained for the 
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full and half electronic spectrum are more important for long chains than for short chains. For the short chains, the 
difference in dimerization amplitude does not exceed « 10%. In the extreme case of the shortest two-atom chain, 
working with half the electronic spectrum gives the exact results. For the long chains, the difference in dimerization 
amplitude increases and is « 17% for the Na = 100 chain length. We attribute the origin of these differences to the 
fact that the gap of the molecular wires decreases with the chain length. However, as we show in the next section, 
an asymptotic regime is reached for chain length Na > 100, where the gap becomes independent of the chain length. 
We therefore can assume that for longer chains, the difference in dimerization amplitude should not increase further. 

We also calculated the corresponding lattice distortions using quantum phonons. The calculations were done by 
determining the ground state of the fully quantum electron/phonon Hamiltonian Eq. (|^) where we introduced a 
cut-off for the number of possible phonon excitations. Each harmonic quantum phonon mode can contain only up 
to n^^^ quanta. Note again, that the calculations were performed with the n, m sums running only over half of the 
electronic spectrum. Then, for the sums running over the originally empty conductance band states, the ground state 
l^'o) corresponds to the situation where an extra electron has been added to the chain. The equivalent situation 
corresponding to removing an electron is obtained by summing the electronic eigenstates only over the originally 
occupied valence band states and considering among these eigenstates of Eq. (|l|) that with the highest eigenvalue. 

In practice, here we choose to calculate the situation corresponding to a chain charged with one extra electron. 
Due to charge-conjugation symmetry, the results for hole injection will be identical. Once the ground state I^Po) 
of Eq. (]l|) is obtained, we can calculate the quantum average (Sq) = (^'o|(5^|^'o) for the mean displacement of the 
phonon mode q, where Sq is given by Eq. ( |A7| ). The atomic displacements induced by charging are obtained 
from (ui) — i^q)- The resulting dimerization patterns are shown on Figure ^. Convergence of the results 

is obtained for a small and finite number of quanta in each mode, roughly n™^'^ w 4,5,6. We will show in the 
next section that the results for the transport properties of the molecular wire coupled to the electrodes converge 
faster with respect to the values of n^^^, especially when one considers charge transport in the tunneling regime. 
As expected, the dimerization patterns are very close to those obtained from the classical phonon model with half 
electronic spectrum. The slight differences may be due to quantum delocalization of the eigenstate of Eq. (^. 

Finally, it can be noticed that the general shape of the dimerization induced by charging and the spatial extent of 
the corresponding polaron arc well reproduced by the different approximations (half electronic spectrum, limited set 
of optical phonon modes, finite number of phonon excitations) introduced to reduce the Hilbert space. The reduction 
of the parameter space is therefore entirely justified. This reduction permits us to treat a great range of molecular 
wire lengths and to determine the transport properties of the wires presented in the next section, with reasonable 
computing times. 



III. RESULTS 



In this section, we present results for the electron transport through typical heterojunctions made of a molecular 
wire connected to two electron reservoirs. In our calculations, the dynamical (i.e. energy-dependent) correlation 
between the electron and the phonon degrees of freedom is kept. Such quantum coherence between the electron and 
the lattice is essential to treat the propagation of polarons through the wire in the tunneling regime {i.e. for an 
injection energy E inside the gap of the molecular wire). We use the term coherence because there is no random 
dephasing of the wave functions introduced by the electron-phonon interaction. Instead each wave vector \j, {uq}) 
has both a definite amplitude and phase, which arc both dependent on the energy £l£3. 



A. Lattice distortions induced by a tunneling electron 

In order to analyse the response of the molecular wire lattice to a tunneling electron in the stationary state, we 
calculate the expectation value of a correlation function between the phonon field displacements and the electron 
density. Such a correlation function is defined as the following quantum average 




where Pi — cjci is the electron wave function projector onto atomic site i (the electron density operator on site 
i). The quantum average of any electron and/or phonon operator O is given by (O) = {0{E)) ~ {"^ {E)\0\i! {E)) . 
The correlation function Jg' represents the mean displacement of the phonon mode q when the electron is on site 
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i. We can then define the conditionally averaged atomic displacement Xj = '^gVq{j) 5q representing the atomic 
displacement on site j on the condition that the electron is on site i. From these atomic displacements, we can obtain 
the dimerization pattern d^*' = (— 1)^ (a;^*|-^ — 2a;^'' + a^jli)- 

We have calculated the response of the molecular lattice to the tunneling electron for different wire lengths and 
different injection energies E. For example, Figure^ shows a two-dimensional map (atomic position j'/electron position 
i) of the dimerization Sj^ for a Na — 100 chain length. The tunneling electron is injected at mid-gap {E — 0.0) from 
the left {j = 1) and propagates to the right (j = 100). The bright part around the first diagonal in Fig. ^ represents a 
dip in the dimerization pattern along the atomic positions. As shown in the previous section, such a lattice distortion 
is the signature of the formation of a polaron inside the chain. We name it a virtual polaron because of the transient 
nature of the tunneling electron injected inside the molecular wire gap. The induced atomic distortions are located 
around the electron position i (vertical axis). As can be seen in Figs. ^ and|^, the virtual polaron has its own intrinsinc 
width (estimated around « 15 atomic sites iar the present model of molecular chains which is consistent with the 
correlation length ^ of the continuum modclEiJ). Its formation is therefore not possible for the electron positions at 
the ends of the molecular wire. For short chains {Na < 15), the lattice distortion associated with the polaron cannot 
be accommodated in the wire. This has an important consequence on the electronic transport properties, as we shall 
see in the next section. 

The atomic displacements are, for most electron positions, slightly less than those for an isolated chain (Fig. ||) 
because the lattice does not respond fully to the tunneling electron as in the case of the static charge added into the 
chain. The width of the lattice distortion around the electron is also smaller for the tunneling electron than for the 
static charge. 

The amplitude of the virtual polaron slightly increases when the injection energy E increases above mid-gap. Such 
a behaviour persists until E reaches the first resonance peak in the electron transmission (see next section) . This may 
be understood from the fact that far from mid-gap, the electron wave function gets more weight (the amplitude of 
the electron wave function gets larger) leading to a stronger coupling to the phonons and therefore to more important 
lattice distortions. Above the first resonance peak, the electron is not, strictly speaking, tunneling anymore through 
the molecular gap. For these energies, the transport regime is better described as tunneling resonantly through the 
coupled e-ph states of the system. The scattering wave function acquires more weight from these quasi-standing 
waves which in turn modify qualitatively the polaron- like nature of the corresponding atomic, distortions. A detailed 
analysis of such distortions is out of the scope of this paper and will be presented elsewhereEil. 

The coherent electron-lattice distortion leads to a modification of the electronic spectrum of the molecular wire 
compared to the spectrum of the undistorted chain. In all the cases studied, the polaron formation is associated with 
a reduction of the gap of the originally undistorted chain. Because the atomic distortions induced in the case of a 
static electron added in the chain are different from those obtained for a tunneling electron, we expect the electron 
transport properties through the corresponding spectra to be different. The behaviour of the transport properties in 
the molecular wires is shown and analysed in terms of electron transmission probabilities in the next section. 



B. Transport properties 

The current flowing through the different channels (for example in the right outgoing channels) is given by 

2e / \ 

jfmj(^) X ^"^ a/+l,{m,}j , (22) 

for sites I > Na + 1 located inside the right lead. From the asymptotic form of the wave functions inside the leads Eqs. 
(^ and (^, it is found that y is related to the transmission probability by j = 2|-/3r sin fc^^^j |f{-„_^i,(i?)p. 

Similarly, the current inside the left lead channels j^^ j is related to the reflection probability |r{-„j_j}(£')p. 
We can define an effective total transmission probability as 

..^ , /3r sin fc? 1 

T{E) ^ hm,m\' a ■ II' ■ (23) 

This takes the form of a sum of contributions from the different outgoing channelsEl. 

In the following, we present results for the injection of an electron and transport by tunneling effect inside the gap 
and by resonant tunneling through the levels of the system above the gap. The electron injection energy is defined as 
positive with respect to the reference energy E = for all the molecular wires. The Fermi energies of the leads are 
assumed to be pinned at mid-gap in absence of any applied bias. 
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Figure ^ shows typical results for the effective total transmission T{E) and the contribution from the elastic channel 
and from one inelastic channel. The different transmission probabilities have, as expected, an exponential behaviour 
in the tunneling regime (below the first resonance peak) and present resonances through the spectrum of the system 
for higher energies. 

It is interesting to note that in the tunneling regime when the injection energies are below the first resonance peak 
{E < 0.51 eV for Na = 40 and E < 0.39 cV for Na = 100), the most important contribution to the total transmission 
T{E) comes from the elastic channel. Above the first resonance, the contribution of the other inelastic channels may 
become as important as the contribution of the elastic channel. Although the contribution of the elastic process is 
dominant in the tunneling regime, it should be noted that this contribution is quite different to that obtained from an 
elastic scattering treatment of the transport through the undistorted molecule. This is because the 'elastic channel' 
corresponds to processes in which, once the electron has crossed the molecule, there is no overall absorption or emission 
of phonons. Nevertheless, phonons are emitted and absorbed during the intermediate stages of the transport, and 
therefore this 'elastic channel' is quite diffeifent from a rigid-molecule calculation where no coupling to the phonons 
is present. We have already shown in RcfJ13 that the transmission is actually enhanced (in the tunneling regime 
and in the limit of low temperatures) because of the e-ph coupling. The transmission enhancement comes from an 
effective gap reduction associated with the virtual polaron formation. This gap reduction is also different from the 
one obtained by charging a (classical lattice) chain with a static electron. 

Figure [t] represents half the HOMO-LUMO gap of isolated neutral and charged chains obtained from the SSH 
model and the effective gap of our coupled quantum electron- phonon model for different chain lengths. Half the gap 
corresponds to the value of the lowest unoccupied molecular states (LUMO) of the isolated (neutral or charged) chain. 
The effective gap of the molecular wire (including e-ph coupling) connected to the leads is obtained from the position 
of the first resonance peak in the transmission T{E). We also plotted the position of the first resonance peak obtained 
from purely elastic calculations (ignoring the e-ph coupling). As expected, the positions of these resonances reproduce 
the values of half the gap of the corresponding neutral chains. An energy shift appears in the first resonance position 
because of the real part of the embedding potentials S'" and that appear in the solution of the scattering problem 
Eq. (0). 

For the fully inelastic scattering calculations, the first resonance in T{E) occurs for injection energy E smaller than 
those obtained from the (solely) elastic calculations. The behaviour characterises the effective gap reduction of the wire 
due to the e-ph interaction. The energy position of the first resonance in T{E) is close to the charging energy E'^^^'^^ 
of the isolated molecular wire. The charging energy £''=^^''8 is defined as the difference between the (selfconsistent) 
ground state total energy of the charged chain E+i [u^] (with the distorted lattice positions corresponding to a static 
polaron) and the ground state total energy of the neutral chain E'ol'^i] (with the undistorted, perfectly dimerized, 
lattice positions u°). Some values of i^'^'^^'g for different molecular wire lengths are given in Table |. Differences 
between the values of E'^^'^'^^ and the first resonance in T{E) do occur, they are due to (i) the systematic energy shift 
introduced by the embedding potentials and also to (ii) the differe nces arising from treating the lattice distortions 
classically or with quantum phonons as pointed out in section [ill A . 

To summarise, there are two distinct physical processes that affect the value of the band-gap of the molecular wire 
and therefore the transport properties through this wire. Firstly, the intrinsic band-gap dimishes with increasing 

wire lengths. The values reach a asymptotic regime for a length Na ^ 100. The asymptotic band-gap value and the 
length above which the asympotic regime is obtained depend of the chemical nature of the molecule, i.e. on the SSH 
parameters used to model the molecular chain (cf. Appendix ^ . Secondly, an effective band-gap reduction occurs 
upon charging the molecular wire with a static charge or with a transient tunneling charge. The reduction is due to 
the the coupling between the charge and the lattice leading to the formation of static or virtual polaron respectively. 
Therefore, the electron transmission through the molecular wire increases because of the e-ph coupling for injection 
energies inside the gap. Although, the trends for the gap reduction are similar for the static and virtual polarons, 
they differ quantitatively over the whole range of molecular lengths studied here. 

Finally, it should be noticed that for short wires {Na < 10) the effective gap obtained from the fully inelastic calcu- 
lations converges towards the gap obtained from elastic calculations. Therefore the values of the electron transmission 
are almost identical for both inelastic and purely elastic calculations. As mentioned in section [II A, although the 
polaron cannot be accomodated in very short wires, in these conditions (short tunneling length and large gap), the 
tunneling process itself is too fast to get significant lattice distortions associated with the charge injection. In order 
to illustrate this point, we give in Table |l| the characteristic times associated with the electron and the phonons. We 
calculate the time domain Tph of the phonons from the extremal phonon frequencies used in our calculations. We 
estimate a residence time t^^s for the electron frarn the full width at half maximum t] of the first resonance peak using 
Tres = h/2r]. Following Biittiker and LandauerOO, a traversal time tbl forlhe tunneling electron can be obtained; it 
is calculated from the elastic channel transmission coefficient as done in Ref.Ej. The traversal time tbl is by definition 
dependent on the electron injection energy E. We give in Table || the range of tbl corresponding to energies inside 
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the tunneling gap. We can see that for short wires {Na <C 2^), the tunneHng time (and also Tres) is smaller than the 
characteristic times associated with the phonons. The lattice dynamics is therefore not fast enough to respond to the 
tunneling electron. In this regime, the transport properties obtained from purely elastic scattering (rigid lattice) are 
not strongly different from those obtained by inelastic scattering (coherently distorted lattice by the e-ph coupling). 
For longer wires {Na ~ 40 2^), the tunneling time (and Trcs) is comparable to or larger than the characteristic 
phonon times; the polaron can be formed inside the molecular wire. It is in this regime that we observe the most 
important differences between the electron transmission for a rigid lattice, and for a lattice that can be deformed by 
the tunneling electron. 



C. Lattice fluctuations 



The imparrtance of lattice fluctuations on the electronic structure of conjugated molecules has already been 
consideredE^lila. As the lattice fluctuations (even in the limit of zero temperature, i.e. the zero-point motion) 
are of the same order of magnitude than the distortions induced by charge injection, it is important to know if such 
species (polarons, solitons) survive the lattice fluctuations. 

We present here results for the electron transmission through a molecular chain where disorder is introduced due to 
lattice fluctuations. In order to compare these results with the full quantum inelastic transmission, we consider the 
limit of low temperatures where the lattice fluctuations are due to the zero-point motion of each phonon mode. The 
averaged transmission is obtained by summing up the different elastic transmission probabilities T{E; { Ag}) associated 
with the displaced phonon conflgurations {Ag}, weighted by the Gaussian distribution probability of the ground state 
harmonic oscillator wave function. The average of a function /({Aq}) depending on the phonon displacements {A^} 
is calculated as 



f exp(-iM[j„A^) 
(/)dis = / dAi...dA,...dA^^„ /({AJ)n ' ' ^ ' , (24) 



where the width of the Gaussian distribution is given by the virial theorem for the zero-point fluctuation of each mode 
iMcj2(x2) = i£;„=o = \{n + \)h.ujq = ^q and Uq = {X^}'/^. 

When /({Aq}) = J2(E] {Aq}), Eq. ( p4[ ) is equivalent to the static lattice approximation expression obtained by 
Pazy and LaikhtmanO using a path-integral formalism. In our calculations, the transmission T{E; {A^}) is calculated 
from elastic scattering through the one-electron spectrum of the molecular chain distorted by the zero-point lattice 



fluctuations. The distorted molecule eigenstates are obtained from the Hamiltonian Eq. (Al) by replacing the 
equilibrium atomic positions by Ui = + Sui where 5ui = '^qVq{i)Aq for the different configurations {A^} of 
the distorted lattice. In practice, T{E; {A,}) is calculated in a similar way as described in section [I A but using 



the approximation of elastic scattering (i.e. no energy exchange between the electron and the phonons is allowed, 
setting n™^^ = for all the phonon modes considered). Futhermore, the integrals over the phonon displacements are 
performed using an algorithm to generate random deviates with a normal Gaussian distribution with a given zero 
mean value and a variance tjq. 

The transmission Eq. (^^ corresponds to an average over the phonon modes, treated classically (as in Appendix 

^) but with a mean-square displacement equal to the quantum zero-point motion, while the electronic transmission 

T{E; {Aq}) is obtained from quantum mechanics. In this average the influence of the phonon displacements on the 
electronic spectrum is taken into account but not vice versa. This is the fundamental difference with the transport 
calculations presented in section [II B| where both the electron and phonon degrees of freedom are treated at the 



quantum level and where the lattice distortion is induced by the injected electron. We therefore expect the results 
for zero-point motion to be different from the one due to the quantum coherent electron-phonon coupling. 

Figure @ shows the elastic transmission probability through the rigid undistorted lattice and the transmission 
averaged (in the manner discussed below) over the zero-point fluctuations. The total effective transmission obtained 
from the inelastic scattering calculations is also shown. 

As expected, the lattice fluctuations substantially modify the electron transmission through the molecular wire. 
The transmission is reduced around the resonance peaks of the one-electron spectrum of the undistorted chain, and 
the peaks are smeared out. However, there is no shift in the position of the resonances corresponding to the formation 
of virtual polarons. 

As in the fully quantum calculations, the transmission is enhanced for injection energies inside the gap of the 
undistorted chain. However, at this point it becomes important to know exactly how the average over the lattice 
fluctuations in Eq. ( p^ is calculated. It is knownEil that the ensemble average of the transmission probability in 
a disordered one-dimensional conductor is a statistically ill-defined quantity, in the sense that it is dominated by a 
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very small number of exceptional configurations. In other words, the mean is in no way representative of the typical 
transmission probability to be expectedJrom a randomly sampled member of the ensemble. This difficulty does not 
arise if logT, rather then T, is averagedB. In Figure ^ we therefore compare logT from our fully quantum results with 
(logT)dis- In this case we see that fully including the electron-phonon coupling (allowing the possibility of polaron 
formation) gives a greater increase in the tunneling transmission than averaging over the disorder introduced by the 
zero-point motion. This shows that only the fully quantum coherent electron-phonon coupling calculations (exact in 
the limit of zero temperature) give the correct physics. 

IV. CONCLUSION 

In this paper, we have presented a method to calculate the inelastic effects on electron transport through one- 
dimensional molecular wires, including realistic electron-phonon coupling. The model for the wires is inspired by the 
Su-Schricffer-Heeger model for trans-polyacetylene. The transport through the quantum electron-phonon system is 
solved by means of a multichannel scattering technique where each channel is associated with probabilities for the 
electron to be reflected or transmitted, given the phonon occupation number configuration. 

The results show that the transport in the one-dimensional molecular wires does not occur as in traditional (three- 
dimensional) semiconducting molecular devices. It does not involve the propagation of free electron- like particles, 
but instead is due to the coherent propagation of "quasi-particles" : an electron (or hole) surrounded by a lattice 
distortion, characteristic of the formation of an electron (or hole) polaron. This object is called here a virtual polaron 
because of the transient nature of the timncling electron (or hole) injected inside the HOMO-LUMO gap of the 
molecular wire. In this regime and in the limit of low temperatures, the tunneling transmission probability through 
the device increases, due to the electron-phonon coupling, in comparison with the transmission obtained from elastic 
scattering through the undistorted molecular wire. Lattice fluctuations also modify the electron transmission through 
the wire. However, the corresponding enhancement of the transmission in the tunneling regime is less than that 
produced by the virtual polarons. 

The influence of other defects (such as a soliton-like defect existing in the chain) on the transmission has already 
been considered, and will be presented elsewhere. The transport for finite temperatures (different initial phonon 
occupations and proper statistical averages) is under study and results will be presented in the near future. 

The results presented in this paper are general and can be applied to other types of one-dimensional atomic-scale 
wires subject to a Peierls transition. For example, it has been shown both experimentally and theoretically that—a 
Peierls-like transition also occurs in dangling-bond (DB) lines fabricated on the H-passivated Si(OOl) surfaceEj E3. 
More recently, it has been shown theoretically that the injectiep of a static charge in the bands around the band-gap 
leads to a distortion of the atomic positions along the lineE^. This distortion corresponds to the formation of a 
small polaron in the DB line. We therefore expect that carrier injection in the band-gap of the DB lines will lead to 
similar physical results to those presented here for molecular wires. Many other quasi-one-dimensional systems may 
be expected to show similar characteristics. 
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APPENDIX A: ISOLATED MOLECULAR WIRES 
1. The SSH model 

Su, Schrieffut-and Heeger (SSH) modelled a chain of trans-polyacetylene (t-PA) as a purely one-dimensional atomic 
(CH)a; chainHH. The model combines a classical ball-and-spring term for the distortion of the u-bond backbone 
with a tight-binding representation for the delocalised 7r-electron orbitals along the chain. Furthermore the electron 
hopping integrals between adjacent (CH) groups are expanded linearly about some reference values. The corresponding 
Hamiltonian is 
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HssH — — 



(to - a {u^+i -Ui)) 



sCi+i,s + a 



Ui) 



(Al) 



where cj^ (ci,s) creates (annihilates) a 7r-electron of spin s at site i. Ui is the displacement of the i-th (CH) group 
from its place in the reference (undimerised) system with hopping integral to- K is the spring constant corresponding 
to the (T-bond and a is the electron-lattice coupling constant. As explained in the introduction, the undimerised 
metallic chain is unstable with respect to a Peierls distortion, and the ground state of an infinite neutral chain has 
displacements given by Ui = (—1)* uq, where uq is a constant depending on to, a and K. The values chosen for the 
different parameters (to=2.5 eV, a=6.1 eV/A and K=A2.Q eV / k^fA give for an infinite perfectly dimerised chain 
a total band width of 10 eV, a band gap of 1.4 eV, and a difference between long/short bond lengths of 0.1 A in 
agreement with experiments. 

For finite size chains containing Na atomic sites (i.e. (CH) groups), the ground state GS(iVa,iVe) of Hamiltonian 
( Al) can be obtained for neutral (A^e = -^a) or charged {Nf, ^ Na) chains, Nf, being the total number of 7r-electrons 



inside thepChain. The ground state is obtained when the restoring forces of the springs balance the electronic forces 
by solvingP^El 



2a 



OCC 

E 



Z„(^) (Z„(^ - 1) - Z„{i + 1)) + K {2u, - u. 



H+l 



)-o. 



(A2) 



where Zn{i) are the compo nen ts of the one-electron eigenstates (with energy e„) of -ffsSH for a given atomic configura- 
tion {ui}. The sum in Eq.(A2) runs only over the occupied states and the spin index s is implicitly taken into account 
in the n summations. For the finite size molecular chains we study here, the boundary conditions are Zn{i) ~ and 
Ui — when i < 1 or i > Na- Furthermore in practice, in order to avoid the uniform translation of the chain through 
space, one atomic site is kept fixed, for instance one end of the molecular chain is kept fixed (i.e. un^ =0). 



2. Classical phonons 

In the limit of small displacements {Sui} around the equilibrium atomic positions {u^} of the ground-. stajtej 
GS(A"a,A"e), it is possible to derive an effective harmonic phonon Hamiltonian by partitioning Eq.( |AlD as followsE3'E3~Ell 

i?sSH[{c„}, {u^ + du,}] = i?.t[{c„}, {u^}] + -ffph[{w°}, {Su,}] + iJe-ph[{c„}, {Su,}] , (A3) 

where Hgt, ^^e-phi ^ph are respectively the static, phonon and electron-phonon coupling parts of the total SSH 
Hamiltonian. _ffe-ph is usually treated as a perturbation up to second order to give pa- quadratic term in Sui in the 
effective phonon Hamiltonian ■ Kij SuiSuj. The dynamical matrix Kij is given bySj 

unocc occ 

K,j = K,, - ^ ^ {F{i, n, m) - F{i + 1, n, m)) {F{j, n, m) - F{j + 1, n, m)) / (e,^ - e„) , (A4) 

n rn 

where F{j,n,m) — Zn{j)Zm{j ^ 1) + Zm{j)Z„{j — 1). The n (m) sum runs over the empty (occupied) electronic 
states. The eigenstates of dynamical matrix Kij give the eigenmodes Vq{i) of vibration (phonons) of the finite size 
molecular chain, while the eigenvalues, Mlu'^, of Kij are related to the phonon frequencies ujq {M being the mass of 
the CH group). 

In the present paper, we use the same boundary conditions as above to determine the phonon modes. However, 
different boundary conditions could be used : fixed ends (i.e. constant molecular chain length), free ends eventually 
coupled to different spring constants to simulate the effective coupling to the electrodes. Although, these different 
boundary conditions would in principle affect the electronic and the vibrational properties of the chairO, we infer 
that the main physical results obtained in the present study will not be drastically modified. For instance, it appears 
that these different conditions will mostly affect the acoustic modes of the chain. Those modes have been however 
neglected in the present work because their contribution to the (virtual) polaron formation is negligible. 



3. Quantum phonons 



At this stage, we already have all the ingredients to derive a quantum version of the SSH Hamiltonian. From 
the reference system, chosen to be the neutral molecular chain of length Na, we can write the Hamiltonians for the 
non-interacting electron and phonon degrees of freedom as 
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n 

where c]j = J^i Zn{i) c\ creates (c„ annihilates) an electron in the nth electronic state of the reference system with 
energy Cn- The harmonic phonon Hamiltonian (neglecting the zero-point energy) is 

where (a^) creates (annihilates) a phonon mode q with frequency usq. The basis set associated to + iJph is 
formed by the eigenstates \n,{nq}) — Yiq ^"y"* i |Q) with eigenvalues = + X]q ^"-^ij' where |0) is the 



vacuum state and {riq} the set of phonon occupation numbers. 

We expand the lattice deformations Sxi induced by an additional charge introduced in the chain onto the phonon 
modes of the neutral chain: Sui — Vq{i) Sq. The new lattice positions, displaced from the equilibrium position u^, 
are Ui = + 6ui. Then, the linear electron-phonon coupling term of the original SSH Hamiltonian is written in a 
quantum form by quantizing the phonon field displacements 



Therefore the e-ph coupling Hamiltonian is 

^e-ph = ^ Iqnm (aj + Og) , (A8) 



q,n,7n 



where 



7qnm = ^ Ag(i) [Z„(i)Z„(i - 1) + Z„(i - l)Z„(i)] , (A9) 



8=2 

and 



Xq(^) = a[Vqi^)-Vq{^-l)]x (AlO) 

The total Hamiltonian for the molecular wire with quantum phonons and linear electron-phonon coupling inspired 
by the SSH model is given by the sum = Hf,-\- iJph -|- ife-ph as in Eq. (|l]) . 

Finally it should be noted that the e-ph coupling matrix elements ^qnm obey some selection rules. Generally, 
Iqnm = uulcss the dircct product Tq ® r„ ® r,„ contains the identity representation (r„, Tm and Tq being the 
irreducible representation of the eigenstate n, m and of Xq respectively). In practice, Z„ and Vq are even/odd 
functions with respect to the center of the molecule. The quantity Zn{i)Zm{i — 1) + Zn{i — l)Zm{i) is even (odd) 
when n -|- m is an even (odd) integer (indexing the eigenvectors Z„ by increasing eigenvalues and Zn=i being even). 
Whenever the quantity under the site i summation in Eq. ( [A9| ) is odd, ^qnm — 0. Although it is not surprising to 
obtain selection rules for the e-ph coupling, their existence is very important in order to reduce the computing time 
of the product i?„|0) needed to solve Eq. (p^). 
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TABLE I. Huang- Rhys factors Sq, the deformation energy AE — -Eof^i] — Eo[Ui], the harmonic distortion energy huiqSq 

and the charging energy E'^^^''^ = i5+i[ui] — -Eoiwi] for different molecular wire lengths. See main text for the definition of the 
different quantities. 
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TABLE IL Characteristic times associated to the electron and phonons. t^cs is the residence time for the electron estimated 
from the width of the first resonance peak, tbl is the Biittiker-Landauer traversal tunneling time for the electron. The values 
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FIG. 1. Schematic representation of the multichannel configurations for a molecular wire connected to two electrodes. The 
different diagrams represent the bond length alternation in the wire for the different channels. Initially the wire is its ground 
state phonon configuration (lower diagram showing a perfect bond length alternation in the middle of the wire) . The incoming 
electron can exchange energy with the phonon modes inside the wire and therefore modify the initial bond length pattern 
(other diagrams). 
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FIG. 2. Lowest energy (longest wave-length) optical phonon modes Vq{j) for the neutral chain containing Na = 100 atomic 
sites. The frequencies (energies) of the modes are hujq=0.136, 0.141, 0.147, 0.153, 0.158 and 0.163 eV for the mode g=41, 44, 
47, 50, 54 and 58 respectively. The phonon mode indexes q of the finite length chain are taken such that the (acoustic and 
optic) phonons are ordered by increasing frequency. 
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FIG. 3. Dimerization pattern dj (in A) induced by adding an extra electron into the molecular wire for three different 
wire lengths A^a = 40, 70 and 100. The dip in the dimerization patterns is characteristic of the formation of a static polaron 
located in the middle of the chain. The dotted lines correspond to the dimerization obtained form the original SSH model with 
the full electronic spectrum. The dimerization obtained by considering half the electronic spectrum and classical phonons is 
given by the solid lines (all the phonon modes), and circles (6 optical modes for Na = 40 and 70 and 10 modes for Na = 100). 
The dimerization calculated for quantum phonons is represented by the dot-dashed lines (6 optical modes and = 6 for 



A^a = 40, 6 modes and n" 



: 4 for Na = 70, 4 optical modes and n^c" = 5 for Na = 100). 



polaron 




atomic position j 



FIG. 4. 2D (atomic position - electron position) map of the dimerization pattern (in A) obtained from the atomic displace- 
100 chain length. For the boundary conditions chosen, an electron is injected from the left (injection 



ments x^^^ for a A'^a 



energy E = 0.0 at mid-gap). The charge current is flowing from atomic site j = 1 to site j = 100 (horizontal axis). The bright 
parts represent a dip in the dimerization pattern, i.e. the formation of a virtual polaron. The induced atomic distortions are 
always located around the (tunneling) electron position i (vertical axis) when there is enough "room" for the polaron to exist 
inside the chain. 
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FIG. 5. Dimerization pattern obtained from the atomic displacements a;''' for Na = 100 (identical to figure ^ for three 
different electron positions i = 19 (solid line), i = 51 (dot-dashed line) and i = 81 (dotted line). The corresponding dimerization 
pattern for a static polaron inside an isolated chain is also shown (thin dashed line, identical to the dimerization shown in figure 
y as the dot-dashed line). Note the difference in the width of the polaron defect for a static electron and a tunneling electron. 
Note also the "wake" of lattice distortion left by the tunneling electron. 
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FIG. 6. Electron transmission probability (on logarithmic scale) versus electron injection energy E through molecular wires 
of different length. Calculations were performed with A'ph = 6, n™^ = 2 for Na = 40 and with Nph = 4, n™^ = 3 for 
Na — 100. Effective total transmission T{E) (solid lines), transmission from the elastic channel |t{o}(£')P only (dotted lines), 
transmission from the first inelastic channel {'n.q=i = 1 and nq>i = 0} (dashed lines). We recall that the first optical mode 
5 = 1 is the optical mode with the longest wave-length {i.e. lowest energy). The corresponding energy is hLOq=0.148 eV (for 
Na = 40) and ?iWq=0.136 eV (for Na = 100). In the tunnehng regime, i.e. for energy E below « 0.51 eV (for Na = 40), and w 
0.39 eV (for Na = 100), the main contribution to T{E) comes from the elastic channel. 
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FIG. 7. Half gap obtained from the original SSH model versus the chain length for a neutral molecule (solid line) and for 

a molecule charged with one extra electron (dashed line). The corresponding charging energy is also shown (dotted line with 
diamonds). The energy positions of the first resonance peak in the transmission are also shown: for purely elastic scattering, 
i.e. ignoring the coupling to the phonons (empty circles), and for inelastic scattering including the coupling of the electron to 
the A^ph = 4 lowest frequency optical modes (with = 2) (filled circles). All energies are given in eV. 
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FIG. 8. Electron transmission T{E) versus electron injection energy E for two molecular wire lengths Na = 40 and 
Na = 100. Transmission through the one-electron eigenstates of the neutral chains from elastic scattering (thin solid lines). 
Transmission through the coupled quantum electron/phonons system from inelastic scattering (solid lines with circles), A^ph = 6 
and n™e'' = 2 for Na = 40, A^ph = 4 and r(™c^ = 3 for Na = 100. Averaged transmission (log r(i5))dis from elastic scattering 
through the one-electron eigenstates of the chain distorted by zero-point lattice fluctuations (dotted lines), Nph = 3 for Na = 40 
and Na = 100. 
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